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M.  ABSTRACT  (CmtHtm  m  wntw  NR  «  MMtwy  m*  tomtit?  h?  MmI  mmmm) 

T  Successful  operation  of  a  tandem  mirror  fusion  reactor  is  critically  dependent  on  electron 
cyclotron  resonance  heating.  We  find  that  for  proof-of-principle  parameters,  this  frequency  is 
more  than  a  factor  of  two  lower  than  the  nonreiathristic  cyclotron  frequency.  Presently 
available  sources  can  provide  relatively  high  power  at  this  frequency,  x 


ELECTRON  CYCLOTRON  RESONANCE  HEATING  OF 
TANDEM  MIRRORS  AT  RELATIVISTIC  ENERGIES 

Electron  cyclotron  resonance  heating  (ECRH)  has  been  shown  to  be  an  efficient  way  of 
heating  plasmas  in  tokamaks  and  Elmo-Bump'y-Torus.1  Recently,  Baldwin  and  Logan2  pro¬ 
posed  selective  local  heating  of  electrons  of  tandem  mirrors,  e.g.,  inside  the  Ying  Yang  and  the 
A-cell  to  lower  the  power  requirement  of  the  neutral  beam  injectors.  Small  scale  experiments 
are  already  being  conducted.3  Wave  propagation  and  energy  deposition  in  tandem  mirrors  are 
also  under  active  theoretical  investigation.4  Most  of  the  research  or  design  so  far  uses  a  heating 
frequency  a*  equals  to  the  first  or  second  harmonics  of  the  cyclotron  frequency  <oc,(0)  at  the 
center  of  the  plasma.  For  a  20  KG  magnetic  field,  this  would  require  a  56  GHz  source,  such  as 
gyrotrons.  High  power  gyrotrons  at  this  frequency  are  still  under  development.  In  the  follow¬ 
ing,  we  show  that  as  the  electron  temperature  ( T, )  increases,  a  lower  frequency  wave  will  actu¬ 
ally  deliver  its  energy  closer  to  the  center.  For  a  20  KG  magnetic  field  and  Tf  -  90  keV,  wave 
frequency  corresponds  to  10  KG  magnetic  field  (28  GHz)  is  sufficient.  Gyrotrons  at  this  fre¬ 
quency  range  are  already  available.5 

To  study  the  propagation  and  energy  deposition  of  the  wave  in  tandem  mirrors,  we  use  a 
ray  tracing  code  which  incorporates  fully  relativistic  absorption  coefficients.  No  expansion  on  y 
is  made,  where  y  is  the  relativistic  factor.  Shkarofsky’s  integrals  are  not  used  here.6  A  detail 
description  of  the  code  and  the  absorption  coefficients  will  be  reported  elsewhere.  Here,  we 
shall  concentrate  on  the  results  of  ECR  heating  of  plasma  in  the  Ying  Yang  region  of  MFTF-B. 
We  shall  also  comment  on  applying  our  results  to  the  A-cell. 
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The  Ying  Yang  magnetic  held  is  an  absolute  minimum  B  configuration  and  is  approxi¬ 
mated  by  the  sum  of  a  mirror  field  and  a  quadrupole  field7 
* 

B  “  Bji/  +  Bp 

where  the  mirror  field  is  given  by  Bw  »  fl-e.  -  1/2  Blri,  and  the  quadrupole  field  by  given  by 

Bp  -  -Kp  B.r(£,  cos  20  -  e»  sin  29) 

where 


l  +  K2(z  +  LJ2)2  1  +  K2(z  -  LJ D2 


and  B.  is  the  derivative  of  B.  with  respect  to  z.  B0  is  the  minimum  magnetic  field  and  is  taken 
to  be  20  KG;  Lm  *  360  cm,  is  mirror  to  mirror  distance.8  The  other  constants  are  Kp  —  0.1, 
a  *  10.74,  Km  —  0.09.  The  density  profile  is  modeled  by 


n(r.z)  -  n0  exp  h(r//r)2  -  (z//z)2} 

where  n0  is  the  density  at  the  center  and  I,  and  /.  are  density  scale  lengths  in  radial  and  z  direc¬ 
tion.  The  values  of  n0,  /,,  /.  are  10n,  48  and  60  respectively.  The  temperature  profile  is 
modeled  by’  Tt(r,z)  -  Te^(n/ n0)‘v'  and  Te0  —  90  keV.  It  should  be  noted  that  the  exact  func¬ 
tional  forms  of  n  and  Te  are  not  crucial  in  our  calculations. 


First  we  use  the  first  harmonic  of  the  extraordinary  mode  with  at  -  ot„(0)  where 
<ur,(0)  —  eBjmc  is  the  cyclotron  frequency  (non-relativistic)  at  the  center  of  the  plasma.  We 
fix  0  —  0°,  r  »  100  cm  and  varies  z  —  0,  50,  100  cm.  We  inject  the  rays  5°  from  the  perpen¬ 
dicular  to  the  magnetic  field.  The  results  are  shown  in  Fig.  1 .  The  rays  are  heavily  damped  at 
the  edge  of  the  plasma.  Reducing  T,  from  90  keV  to  10  keV  does  not  improve  the  picture  too 
much.  If  we  decrease  o>  from  id„.(0)  to  some  fraction  of  otf,(0),  we  find  the  energy  deposition 
regions  shift  toward  the  center  of  the  plasma  (Fig.  2).  Further  reduction  of  at  encounters  a 
mode  conversion  and  the  ray  is  reflected  back.  Contrary  to  the  accessibility  condition  w  >  ai„, 
which  requires  a  higher  at  to  reach  higher  density,  we  find  that  lowering  the  frequency  will  actu¬ 
ally  help  to  deliver  the  energy  to  the  center  of  the  plasma.  Similar  results  are  obtained  for  the 


i 
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ordinary  mode  (Fig.  3),  except  for  this  case,  we  can  only  reduce  u  to  0.5  0*^(0).  Further 
reduction  in  «  results  in  violation  of  the  accessibility  condition  >  <0^  and  the  ray  is  turned 
back. 


The  physics  behind  could  be  easily  understood  from  examining  the  relativistic  cyclotron 
resonance  condition 

<d  -  mjy  +  A„vn. 

To  simplify  the  arguments,  we  shall  approximate  k»  =  0  for  almost  perpendicular  injection. 
The  resonance  condition  becomes  w  =  o>ce/y.  For  a  heating  frequency  a>  -  wtr(0)  and  rela¬ 
tivistic  temperature  (y  >  1),  cott,  must  be  bigger  than  w<(P(0)  to  satisfy  the  resonance  condition. 
The  Ying  Yang  region  is  an  absolute  minimum  B  configuration.  If  the  resonance  condition  is  at 
u>ce  >  the  resonance  will  be  located  away  from  the  center  of  plasma,  for  Tc( 0)  -  90 

keV,  it  is  near  the  edge  as  demonstrated  in  Fig.  1. 

In  order  to  move  the  resonance  closer  to  the  center  of  the  plasma,  we  have  to  lower  the 
magnitude  of  01  at  resonance.  This  can  be  accomplished  by  lowering  the  heating  frequency  w 
to  e*»(.,(0)  where  c  <  1.0.  Fig.  2  shows  <  could  be  as  small  as  0.4.  From  the  resonance  con¬ 
dition,  we  find  that  y  must  be  about  2.  For  Te( 0)  -*  90  keV,  this  means  the  energy  is  depo¬ 
sited  at  the  tail  of  a  maxwellian  distribution  function.  We  have  not  optimized  the  injection 
geometry,  it  may  be  for  some  particular  parameters,  y  could  be  smaller.  A  detail  analysis  of 
the  Fokker  Planck  equation  and  the  quasilinear  diffusion  equation  will  be  needed  to  find  how 
fast  the  energy  is  deposited  at  the  tail  and  how  fast  it  is  transferred  to  the  thermal  electrons. 

The  preceding  scheme  could  also  be  applied  to  the  heating  of  the  A-cell  of  MFTF-B,  but 
it  turns  out  that  lowering  the  heating  frequency  may  not  be  necessary.  The  A-cell  is  not  a 
minimum  B  configuration  in  the  radial  direction,  so  the  wave  with  w  *  w,r(0)  should  be  able  to 
deliver  its  energy  to  the  thermal  electrons  and  close  to  the  center  of  the  plasma.  The  central 
magnetic  field  strength  in  the  A-cell  is  about  10  KG,  so  again  a  28  GHz  source  is  sufficient. 
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In  conclusion,  we  find  that  in  the  relativistic  regime,  the  wave  energy  could  be  delivered 
closer  to  the  center  of  the  Ying  Yang  plasma  of  MFTF-B  if  u  <  <oce(0).  A  scenerio  to  heat  the 
plasma  in  this  region  could  be  the  following:  a  series  of  source  at  28  GHz  and  56  GHz  are  used, 
first,  the  56  GHz  source  is  turned  on  and  then  as  the  temperature  of  the  electrons  rises  above  a 
certain  level,  we  switch  off  the  56  GHz  source  and  turn  on  the  28  GHz  source  to  continue  the 
heating. 

We  wish  to  thank  W.  Manheimer,  E.  Ott  and  J.  Shearer  for  discussions.  This  work  was 
supported  by  the  Department  of  Energy. 
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Fig.  1  -  A  r-2  plot  of  the  Ying  Yang  region  of  MFTF-B.  The  extraordinary  modes  are  injected  at  r  -  100  cm,  0  -  0* 
and  z  m  0*  50,  100  cm,  7,(0)  —  90  keV,  no  “  10n  cm  ~3.B0  “  *0  KG,  K0  -  «/«<,(0)  »  1.0.  The  maiority  of  the 
energy  (>95%)  are  deposited  at  the  broken  curve  regions  (■•*).  The  wriggles  (•**-*)  is  energy  desposition  regions  when 
7,(0)  -  10  keV. 


Fig.  2  —  Same  parameters  as  Fig.  1.  The  broken  curve  (-),  open  circles  (OOOO)  and  wriggles  (w) 
are  energy  deposition  regions  of  Y0  -  1.0, 0.75  and  0.4  respectively,  7,(0)  -  90  keV. 
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